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Abstract 

In the recent papers, we introduced a method utilised to measure the flow field. The method is based on the tracking of 
supergranular structures. We did not precisely know, whether its results represent the flow field in the photosphere or in some 
sub-photospheric layers. In this paper, in combination with helioseismic data, we are able to estimate the depths in the solar 
convection envelope, where the detected large-scale flow field is well represented by the surface measurements. We got a clear 
answer to question what kind of structures we track in full-disc Dopplergrams. It seems that in the quiet Sun regions the 
supergranular structures are tracked, while in the regions with the magnetic field the structures of the magnetic field are dominant. 
This observation seems obvious, because the nature of Doppler structures is different in the magnetic regions and in the quiet 
Sun. We show that the large-scale flow detected by our method represents the motion of plasma in layers down to ~10 Mm. 
The supergranules may therefore be treated as the objects carried by the underlying large-scale velocity field. 
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1. Introduction 

The supergranulation is a convection-like pattern that 

p54| ). 



remains a puzzle since its discovery by Hart 



There are many papers published examining the para 
metric properties of supergranular cells. Typical values 
of the sizes are around 20-30 Mm, lifetime of 24 hours 
(see e.g. iDeRosa and Toomre 2004 and references 
herein). There is no explanation yet to give a clear state- 
ment about the origin of supergranules. The classical 
mechanism invoked to explain the origin of the super- 
granulation is the latent heat of the recombination of 
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He^^ into He^ at roughly 10 Mm below the photosphere 
( |Gierasch| |1985| ). In a fluid at rest, such a heat release 
may trigger an instability, which can turn into motions 
at scales comparable to the supergranular scale. The is- 
sue in the Sun is that the plasma in the convection zone 
is highly turbulent and the mentioned thermal instability 
can be suppressed by turbulent motions. There is also 
a lot of works doubting the existence of the supergran- 
ulation as the convection mode. The surface properties 
of the supergranules may be explained as non-linear in- 
teraction between granules triggered by exploding gran- 



ules (e.g. Rieutord et al. 2000 ). The wave-like prop- 
erties of the supergranulation were also discussed (e.g. 
IGizon etal.l 12003] ). 



The time-distance local helioseismology ( jPuvall 
|et al.[ |1993| ) is a tool that allows to examine the struc- 
ture of the plasma flows beneath the solar surface. The 
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/7-modes of the solar oscillations are excited in the up- 
per convection zone and travel through the convective 
envelope. Plasma motions and variations of the sound- 
speed change the speed of the wave-packet propagation. 
From the difference between the theoretical travel-time 
from one point in the solar photosphere to another one 
and the theoretically calculated one, one can map the 
disturbances causing the travel-time deviations in the 
subphotospherical layers. Initial results of the time- 



distance helioseismology from SOHO/MDI ( Duvall layer. Hathaway et al 



et al. 1997) showed that a pattern of the surface flows 
can persist 2-3 Mm below the surface, improved re- 
sults gave 8 Mm ( |Duvall[[l998| ). Further time-distance 
analyses suggested the depth of the supergranular flow 
to be ~ 15 Mm (Zhao and Kosovichev| |2003) . In that 
study, the authors calculated the correlation of the hori- 
zontal divergence at different depths with the horizontal 
divergence just below the surface. They found a posi- 
tive correlation in the depths of 0-5 Mm. Deeper down 
it changes the sign. It may mean that the return flow 
in the supergranules starts to dominate the intercellular 
flow structure. 

The supergranules are well visible on the entire solar 
disc in Dopplergrams. This is caused by the prevailing 
horizontality of the internal velocity field within the su- 
pergranular cells. If we assume that supergranules are 
objects carried by the flow field on the larger scale, the 
visibility of supergranules in the Dopplergrams gives an 
opportunity to use them to map this underlying large- 



termined from MDI Dopplergrams is 5 % larger than 
the corresponding spectroscopic rate. Some explana- 
tions of this phenomenon appeared in the following 
papers. Gizon et al. ( 2003| ) suggested the modulation 
of the convective pattern by travelling waves. Already 
Wolff| ( p^95] ) noted that the effective viscous dissipation 
of Rossby- waves (r-modes), which happen in the Sun 
at 0.932 Rq (~50 Mm below the photosphere), could 
modulate the convection in the near- subphotospherical 
( |2006l ). 



scale flow field. It was first done by Snodgrass and Ul 



rich ( |1990 ). Recently, we developed a method based on 
the local correlation tracking algorithm (LCT; Novem- 
ber! [1986) with a similar util isation. This method (for 
details see Svanda et al. 2006 1 uses the pattern of the su- 
pergranulation in the SOHO/MDI Dopplergrams ( Scher- 



rer et aH|1995) to measure the motions on larger scales. 
The issue of the method is that we measure displace- 
ments of the supergranular structures in the series of 
the processed Dopplergrams and interpret them as the 
large-scale velocity field in the solar photosphere. How- 
ever, we cannot establish the range of depths in the solar 
convection zone, where the large-scale flows are well 
represented by the surface measurements. With the use 
of the data provided by the time-distance helioseismol- 
ogy we can determine this unknown parameter. We can 
also verify the assumption of our analysis: whether the 
supergranules are subjected to the transport by the ve- 
locity field on the larger scales. 

There is a list of works suggesting that the super- 
granules do not represent the actual differential rota- 
tion or meridional circulation. B eck and Schou| ( [2000| ) 
mentioned that rotation rate of the supergranules de- 



however, used synthetic 
data to explain the apparent super-rotation of the super- 
granulation as the projection effects on the line-of-sight 



Doppler velocity signal. [Meunier and Roudier] ( |2007 ) 
re-analysed Hathaway 's idea and found that the projec- 
tion effects can explain the most of the supergranular 
super-rotation, but not all of it. IMeunier and Roudier 



( 2007| ) concluded that the supergranules are unreliable 
in the analy sis of the motions o n supergranular scales. 
However, in| Svanda et al. ( 2007t the results showed that 
the large-scale apparent motion of supergranules mea- 
sured by the method introduced in Svanda et al. (2006] ), 
and used also in this study, is in almost perfect agree- 
ment with the near-surface mass flow measured by the 
time-distance helioseismology. It is perhaps a k-oj fil- 
ter used for the suppression of the noise, which makes 
this method more useful than the previous ones based 
on similar principles. 

The large-scale flows are a combination of many 
types of motions on various spatial scales (such as the 
differential rotation, meridional circulation, possibly the 
giant cell convection, and others), which unfortunately 
cannot be reliably separated in components. These com- 
ponents may vary independently with time and depth. 
This implies that, in principle, when dealing with the 
large-scale flows, we cannot obtain unambiguous results 
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Fig. 1. A sample of acoustic ray paths used for the time-distance 
helioseismology, shown in a vertical plane. The shadowed regions 
illustrate ranges of averaging. The vertical and horizontal lines show 
a grid used for inversion of acoustic travel time data. After |Koso-| 
[vichev and Duvall| p006; , who used the same dataset as in this 
study. 
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that do not allow any alternative interpretation. Nonethe- 
less, the realistic numerical simulation should reveal the 
combined large-scale plasma flow, properties of which 
should be comparable with the measurements. Such nu- 
merical simulation, which is not present at the time, 
shall allow to distinguish various components of the 
detected large-scale flow and to study them separately. 
The results in this paper therefore provide encouraging 
large-scale flow properties, which are to be reproduced 
by models. 



2. Data and Method 

For the comparison between the data obtained by 
our method and the results of the time-distance helio- 
seismology we used high-cadence full-resolution full- 
disc Dopplergrams recorded in March and April 2001 
by SOHO/MDI. In this period, 46 flow maps of the 
vicinity of the active region NOAA 9393/9433 averaged 
over 8 hours were computed by the time-distance helio- 



seismology using the ray approximation (see Zhao and 
[Kosovichev 2004). From these maps we used five ex- 
amples in five diff'erent days (March 3rd, March 28th, 
29th, 30th, and April 25th). In these days, MDI data 
series did not contain many gaps, so the data were suit- 
able for our tracking method, and the field-of-view was 
far from the solar limb, so we can exclude any possible 
limb eff'ect. 

The one-day series of full-disc Dopplergrams were 
processed fol l owing the procedure introduced by 



Svanda et al. (2006) with a few modifications. First, 



the Postel projection was used instead of the Sanson- 
Flamsteed one, because helioseismic data are computed 
in this geometrical projection (it conserves the main 
circles, therefore is suita ble for the measure ment of 
/7-modes travel-times). In Svanda et al. (2006) we per- 



formed the computations in two main steps. The first 
step was used for the determination of the diff'eren- 
tial rotation profile, which was then removed during 
the remapping of the data to decrease the maximum 
displacement caused by motions of supergranular struc- 
tures. In the second step, the flow field with respect 
to the removed mean rotation profile was calculated 
using the algorithm with a higher precision. Then the 
diff'erential rotation profile calculated in the first step 
was added to the flow field calculated in the second 
step. In this study the computations are done in one 
step, i.e., no diff'erential rotation is removed before the 
application of the finer LCT algorithm. We assume that 
in a small field-of-view, the change in the diflferential 
rotation profile is not important. 



In summary, the data processing routine consists of 
the ;?-modes removal (using the weighted average over 
30 minutes), coordinate alignment (all frames record 
the same region on the Sun), coordinate transformation 
(using the Postel projection) and a k-co filtering with 
the cut-off' of 1500 ms"^ for the noise removal. Then, 
the LCT algorithm for the determination of supergranu- 
lar structures displacements is used to calculate the 24- 
hours averaged large-scale velocity field with the 60'' 
correlation window. 

( |2QQ6| ), 



As we discussed in Svanda et al. 



provides the measurements with a random-error of 
15 ms"^ It is believed that the time-distance helio- 
seismology errors are within 10 % in the near- surface 
layers and generally increase with depth, under some 
20 Mm the details of the results must be considered un- 
reliable (|Zhaoj[2006 ). Unfortunately, it is impossible to 
determine the actual errors using the ray approximation. 
The study made by [Birch and Felder| ( |2004| ) showed 
that the estimated errors in the supergranular or moat 
ffow achieved by both ray and Born approximations 
are of the order of 10 %. The development of the helio- 
seismic method is fast and we can expect more reliable 
results with exactly determined accuracies using a per- 
turbation theory and fully consistent OLA inversions 



( Jackiewicz et al.|[2008| ) in a few years. 

All the datasets were aligned with the centre of the 
field-of-view on Carrington coordinates / = 148.5° and 
b=\9\ The field-of-view has the total size of 512x512 
pixels with the resolution of 0.92'' px"^. The datacubes 
from the local helioseismology contain 15 irregularly 
spaced depths from 0.77 Mm to more than 80 Mm. The 
vertical resolution is based on the ray approximation. 
The numbers rather than isolated levels represent the 
consecutive intervals in the radial direction, in which 
the results are averaged. Fig. [T] shows examples of the 
acoustic ray paths, which clearly illustrates the vertical 
resolution. 

The results of our tracking method have an eff'ective 
resolution of 60", but the results of the time-distance he- 
lioseismology have 8" . We binned the helioseismic ffow 
maps to match the resolution of 60". This also means 
that the structure of the internal ffow in the supergran- 
ulation is ffltered out and does not disturb the study of 
the large-scale ffow field. This is the most important de- 



viation from the study of |Zhao and Kosovichev] ( [2003] ). 

The velocity field obtained by our tracking method 
is in principle twodimensional - horizontal. The z- 
component of the time-distance ffows suff'ers from ef- 
fects of cross-talks, therefore we use only horizontal {x 
and y) components. For each depth in the helioseismic 
datacube we calculated its similarity to the ffow map 
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Fig. 2. The dependence of the similarity of the velocity fields provided by the time-distance helioseismology and the flows calculated by the 
LCT method on the depth in the solar convection zone. Left: in quiet Sun regions, right in the regions occupied by the magnetic field. 



obtained by the tracking of supergranular structures. 
As the measure of the similarity we used the magni- 
tude weighted cosin e of the direction difference (as in 
Svanda et al. 2007 ), which is robust to the presence of 
the noise. This quantity is given by the formula: 



Zl^l 



d-b 



Pw 



Zl^l 



(1) 



where a and b are vector fields, a-b is 3. scalar multipli- 
cation and \a\ is a magnitude. The closer this quantity is 
to 1, the better the alignment between two vector fields. 
In our case, a = vlct and b = Vtime-distance (depth). The 
magnitude weighted cosine is also stable to any possi- 
ble issue involving the magnitude determination error 
in the time-distance inversion. 

We investigate the similarity of both methods results 
separately for the regions occupied by the magnetic field 
and the regions of the quiet Sun. Two regimes (magnetic 
and non-magnetic) are separated in the field-of-view 
applying the masks calculated from the surface MDI 
magnetograms. We smoothed the magnetograms by 20'' 
and estimated the threshold of 5 Gauss in the smoothed 
magnetograms as the level distinguishing the magnetic 
regions from the non-magnetic ones. 

3. Results 

3.1. Large-scale horizontal flows 

The results of the comparison are displayed in Fig. [2] 
We basically see a large similarity of the velocity fields 
measured by both methods at higher levels of the solar 
convection zone. The maximum is reached in all but one 



studied cases at the level of 4.5-6.4 Mm (~ 0.006 Rq). 
For the quiet Sun regions, the absolute value of the 
magnitude weighted cosine is not significantly diff'erent 
in the layer 4.5-6.4 Mm and in the neighbouring layer 
3.0-4.5 Mm, it is within the estimated 10% accuracy. 

Our analysis shows that the large-scale flows do not 
vary much with depth down to 10 Mm in the quiet Sun 
regions. This fact means that the large-scale horizontal 
velocity field measured on the surface represents well 
the large-scale horizontal dynamics in the layers of 0- 
10 Mm in the depth. See Fig.|3]for examples of the flow 
maps at difl'erent depths. We cannot expect much better 
match, because both methods measure difl'erent proper- 
ties. While LCT measures the motion of structures, the 
time-distance helioseismology maps the plasma flow. 
Both results may be very difl'erent if any wave-like phe- 
nomenon is involved. 

In magnetised regions, the dominant structures in the 
Dopplergrams are the structures connected rather to the 
magnetic field than to the supergranulation. This idea is 
supported by our results, where we see that in the mag- 
netised regions the detected flows are coherent within 
the depths down to ~5 Mm with the peak in the depths 
of 4.5-6.4 Mm. This is the depth, where downward con- 



vergent flows are still detected is sunspots (e.g. Zhao 



et al. 2001 ) and where sunspots stop exist as the re- 
gions of the suppressed heat transport. The existence of 
the Evershed flow in sunspots, which is easily seen in 
Dopplergrams, supports our interpretation of the veloc- 
ity field in magnetised regions, too. 

There is an alternative interpretation of the obtained 
results. Assuming that the time-distance helioseismol- 
ogy is less reliable deeper down in the convection zone 
and that the large-scale flows do not vary much in the 
upper convection zone, our results could put the level. 
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Fig. 3. The examples of the flow maps in the vicinity of active region NOAA 9393 on March 29th 2001. We see that the topology of flows 
measured by the local helioseismology is similar with the surface one measured by LCT down to depth of roughly 10 Mm, below that the 
topology changes. The lengths of arrows in the time-distance flow maps are not in the same scale, their magnitude increases with depth. In 
the "Supergranular motion" subframe the contours of the mask distinguishing the magnetic and non-magnetic regions are overplotted. 



where the rehabihty of the time-distance measurements 
is reasonable. In this interpretation, the time-distance 
flow measurements are reasonable in upper 10 Mm in 
the quiet Sun regions and upper 5 Mm in regions oc- 
cupied by the magnetic field. There are recent papers 



(e.g. |Jackiewicz et aL||2Q08| ) using another time-distance 
implementation, which suggest that the current time- 
distance results are reasonable only in the upper-most 
3 Mm. 



3.2. Vertical structure of the super granulation 

The obtained results can be interpreted as a fulfil- 
ment of the basic assumption of our tracking method: 
The supergranular structures may be indeed treated as 
the objects carried by the large-scale velocity field. The 
supergranules move in the depths down to 5 Mm coher- 



ently in the horizontal direction with a slight loss of the 
coherence between 5 and 10 Mm in the depth. To con- 
firm this interpretation, the structure of the supergran- 
ules itself must be analysed. 

Therefore we modified our analysis and instead of 
the smearing the time-distance flow field to remove the 
signal of the supergranular internal mass flow, we sub- 
tracted the large-scale velocity field, so that only the in- 
ternal velocities in supergranules remained in the flow 
map. From the maps containing only the separated su- 
pergranulation signal we calculated the horizontal di- 
vergence (see Fig. [4]) as the representation of down/up- 
flows. It is more stable than the measured z-component 
of the velocity, which sufl'ers of plenty of issues ( |Zhao| 
[2006] ). 

From the sequence in Fig. |4] we see that the diver- 
gence signal is very similar in the quiet Sun regions in 
the depths 0.77-6.4 Mm. It is confirmed when plotting 
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Fig. 4. Horizontal divergence maps of the internal supergranular flow at various depths in the solar convection zone. Dark shades of grey 
denote inflows, while light shades represent outflows, shades of grey are scaled automatically for an easy comparison of the structure between 
each other. 



the correlation of the near- surface divergence map with 
maps at various depths (see Fig. |5]). Down to ~7 Mm 
the correlation is positive, implying the coherence in the 
structure of the supergranules. Deeper down the corre- 
lation turns negative, which might suggest the evidence 
of the return flow in supergranules. The typical spatial 
scale of divergence structures goes down by the factor 
of 2-3 in layers deeper than ~8 Mm. A very similar 
analysis was done by |Zhao and Kosovichev| ( |2QQ3| ) on 
another dataset with not much diff'erent results. 

From our data, we cannot compare the structure of 
the supergranulation between the surface measurement 
(the Dopplergram) and the top-most helioseismic layer 
(0.77-1.77 Mm). This topic was, however, assessed by 
Gizon et 3T]pQ00| ). Authors of that paper demonstrated 
the very high correlation between the measured surface 
Dopplergram and the reconstructed line-of- sight veloc- 
ity image from the time-distance flows in the depths of 
0-2 Mm. Based on that and this paper we may conclude 
that the supergranules are highly coherent in the depths 
of 0-6.4 Mm. 

In Fig. |5] the trend of the superadiabaticity within the 
sub-surface layers of the convection zone coming from 
the reference solar model S (|Christensen-Dalsgaard 



et al. 1996) is overplotted. The superadiabaticity A* is 



defined as 



A\r) 



1 d In p(r) d In p(r) 



(2) 



y(r) d In r d In r 
where y, p, and p are the state parametres (adiabatic ex- 



ponent, pressure, and density) of the plasma at distance 
r from the centre of the Sun. The layers are convectively 
unstable, where A*(r) is negative (see discussion in |De-| 
marque et al.[ |1997| ). It is to be noticed that A* turns 



more negative at the depths of 10-12 Mm, where it 
is assumed that the supergranulation should be formed, 
and remain very negative up to the solar photosphere. 
It is the same layer, where the large-scale flow is co- 
herent. This observation may be interpreted that the su- 
pergranules, which form at the depth some 10 Mm be- 
low the surface, are carried by the large-scale velocity 
field, which operates here (and perhaps deeper). There- 
fore the large-scale flows remain nearly constant within 
upper 10 Mm of the convection zone and are detected 
using our tracking method in surface measurements. 



4. Conclusions 

Let's put two separate results obtained in the previ- 
ous Section together. Using the comparison of the large- 
scale horizontal flows obtained by our method based on 
the tracking of the supergranular structures and by the 
local helioseismology we state that the results of the 
tracking method represent the horizontal dynamical be- 
haviour in the depths down to -10 Mm below the solar 
surface. We believe that this result may be important for 
those, who develop numerical models of the large-scale 
convection. 

The supergranulation is coherent within layers, which 
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Fig. 5. Correlation of divergence maps from Fig.|4]at various depths 
with the surface layer in regions of quiet Sun. With the thick line, 
the superadiabaticity in the layers calculated from the model S 
jChristensen-Dalsgaard et al^|1996t is overplotted. 



show correlated large-scale horizontal flows. The inter- 
pretation of supergranules is closer to the real convec- 
tive cells carried by the underlying velocity field. The 
depth of the loss of coherence in large-scale horizontal 
flows coincides with the layer, where the supergranules 
as a strong pattern should be formed. 

The results on the vertical structure of the supergran- 
ulation we show in this paper, however, cannot be pre- 
sented as the general results for the supergranulation 
all over the Sun. We have to keep in mind that these 
results were obtained only using one suitable dataset 
in the quiet Sun regions near the very large active re- 
gion. Much more datasets are needed to describe the 
behaviour of the supergranulation at difl'erent stages of 
the magnetic activity. Perhaps the helioseismic data ex- 
pected from the Solar Dynamics Observatory (SDO) 
will provide some new clues in this topic. 

The ray approximation used in the computation of the 
helioseismological data is often in doubts, suggesting 
that better (e.g. Bom approximation) approach is nec- 
essary. It is known (e.g. |Birch et aL) |2Q01| ) that the ray 
approximation underestimates the small-scale perturba- 
tions. The real perturbations inside the Sun are likely 
to be stronger than those inferred using the ray theory. 
Inversions, however, involve complicated averages of 
observed travel times and are regularised. As the result, 
the details of the inaccuracies of the ray approximation 
on ray-based inversions are not yet clear. Some com- 
parative studies (e.g. [Jensen et al.|[2001| ) comparing the 
results of the ray and Born approximation approaches 
showed that for the large-scale flow they reasonably 
match. It was argued that the ray approximation could 
give credible results with fewer computations. The struc- 
ture of the flow on supergranular scales was tested in 
fZhao et al. (2001 ) with satisfying results. 

At this point we cannot exclude that our results in- 



stead of bringing the new information in the behaviour 
of the supergranules rather put limits on the reliability 
of the time-distance helioseismology. If this interpreta- 
tion applies, then the time-distance measurements of the 
large-scale flows are reliable within upper 10 Mm in the 
quiet Sun regions and within upper 5 Mm in the magne- 
tised areas. Only the improvement of the time-distance 
helioseismology methods can resolve this ambiguity. 
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